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Abstract 

In order to assess the ability of jet observables to constrain the characteristics of the medium produced in heavy-ion 
collisions at the LHC, we investigate the influence of background subtraction and jet quenching on jet reconstruction, 
with focus on the dijet asymmetry as currently studied by ATLAS and CMS. Using a toy model, we examine the 
influence of different background subtraction methods on dijet momentum imbalance and azimuthal distributions. We 
compare the usual jet-area based background subtraction technique and a variant of the noise-pedestal subtraction 
method used by CMS. The purpose of this work is to understand what are the differences between the two techniques, 
given the same event configuration. We analyze the influence of the quenching effect using the Q-PYTHIA Monte 
Carlo on the previous observables and to what extent Q-PYTHIA is able to reproduce the CMS data for the average 
missing transverse momentum that seems to indicate the presence of large angle emission of soft particles. 
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1. Introduction 

One of the main subjects in heavy-ion collisions at ultra-relativistic high energy is the characterization of the 
produced medium through hard probes, like the suppression of high momentum particles or jet-like correlations. With 
the advent of the LHC, jet observables |[T] [2] have also moved to the focus on the discussion. In order to extract 
precisely the medium properties, the analysis of jet production requires several components: generation of medium- 
modified jet events together with a realistic background; and reconstruction of jets and background subtraction as 
close as possible to the experimental analysis. But this program faces several difficulties, like how to model and 
generate a realistic fluctuating background and how to treat the coupling between the hard event and the background. 
Additionally, one needs to understand the different ingredients used in a jet analysis and how they affect the different 
observables. In this work, we want to investigate the effect of background subtraction and of quenching on several 
jet observables. For that, we use Q-PYTHIA[4| jets embedded in a simulated background, which follows a thermal 
plus power law p ± spectrum. In order to easily tune the background parameters, we choose to use a toy model rather 
than a more realistic one. This paper is structured as follows: in section [2j a summary of the background model, jet 
reconstruction and jet subtraction methods is provided. In section|3]a comparison between the different jet subtraction 
techniques for different jet observables is shown, followed by the effect of a jet quenching model, as implemented in 
the Q-PYTHIA Monte Carlo in section [4] Finally, the conclusions are presented in section [5] 
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2. Background Model and Jet Reconstruction 

The background model that is used in this study is based on a thermal spectrum, where the inverse slope, T, the 
only free parameter of this model controls the exponential decay of soft particles. The hard part of the spectrum is 
constrained by a power law, p T b . Particles are generated uniformly in (77, <p) and the multiplicity per pseudorapidity 
unit is fixed to 2100 ± 5% particles. The temperature is chosen so that the region-to-region background fluctuations 
are of similar level as the ones measures by ALICE|5| (corrected to include neutral particles), although the average 
level of contamination is somehow larger. Hence, the chosen values were T = 0.7, T = 0.9 and T = 1.2 GeV, that 
correspond to an average level of fluctuations of 7.69, 10.75 and 15.14, respectively. For jet reconstruction we use 
the anti-kt algorithm with a radius R = 0.3. For the jet background subtraction we compare the standard FastJet[6| 
method based on jet areas using the kt-algorithm with an R — 0.4 over a full stripe in \rj\ < 2, and a variant of 
the CMS technique[7], which is based on a "noise/pedestal" subtraction. In this case, the phase space is divided 
according to the CMS calorimeter segmentation. Then, an average transverse energy, < Ej Ke '\rj) > and a dispersion 
a J j wer (ri) is computed for each stripe in 77. To each cell, the corresponding < E t £ v>er (rj) > and <Tj wer (T]) is subtracted 
and set to zero in case of a negative value. Using only particles that are in positive cells, jets are reconstructed with 
a minimum transverse energy above a cut, Ejj ets . All towers that are contained inside this list of jets are excluded 
and new background parameters are computed with the remaining list of particles. These new parameters are the final 
background estimation values used to perform a final correction to the initial transverse energy of each cell. Again, 
using only non-zero towers, the final list of jets is found. Since the final estimation of parameters should be the close 
as possible of the true background parameters, the parameter Erj e ts> that according to CMS is 10 GeV, changes with 
the temperature, in the case of our toy model. In the end, the best estimation of the background parameters is given 
by E Tjets = 40 for T = 0.7, E Tjets = 60 for T = 0.9 and E Tjets = 70 GeV for T = 1.2 GeV. 



3. Jet Subtraction on Jet Observables 
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(a) Dijet Asymmetry using jet area method (b) Dijet Azimuthal Correlation using jet area (c) Dijet Azimuthal Correlation using CMS- 
for background subtraction. method for background subtraction. like method for background subtraction. 

Figure 1 : Jet Observables for a simulation using Q-PYTHIA with a q = GeV 2 fm~' embedded in a background with different temperature 
parameters. The red/dotted curve corresponds to a T = 0.7, the black/dashed to a T = 0.9 and the green/solid to a T = i.2. The blue points are the 
CMS data and the purple triangles the CMS Monte Carlo simulation. 



In this section, our aim is to compare both background subtraction techniques that were outlined in the previous 
section, and to understand what is the effect of the background parameters, specifically the fluctuations, on the ob- 
servables that we want to describe, dijet asymmetry and dijet azimuthal correlation. Here we present a comparison 
between the two background subtraction methods: CMS-like and FastJet (jet areas). In figure 1(a) it is shown the 
effect of the fluctuations on the dijet asymmetry when using a background subtraction based on jet areas and as one 
can observe the effect of the fluctuations goes in the same direction than data, but no "realistic" values can account 



L. Apolindrio, N. Armesto and L. Cunqueiro /Nuclear Physics A 00 (2012) i-^] 



3 



for the large asymmetry indicated by CMS data. On the contrary, the angular reconstruction of the jets seems to be 
unaltered, as one can see from figure 1(b) As for the CMS-//A:e method, by construction, the resulting dijet asymme- 
try has a much smaller sensitivity to fluctuations, since the jet energy reconstruction is more accurate, but presents a 
larger deviation for the azimuthal correlation (figure 1(c) I. 



4. Quenching on Jet Observables 

Now we turn to evaluating the degree of quenching of the data using Q-PYTHIA. For that, we investigate the 
effect of the influence of the q parameter of this model, by changing its value up to 8 GeV 2 far 1 while keeping the 
background fixed to cr — 10.75 GeV. By increasing q in Q-PYTHIA, the jet shower is accelerated with respect to 
vacuum due to the increasing number of medium-induced splittings. In figure [2] one can observe the effect in the 
dijet asymmetry and azimuthal correlation using jet area method to subtract the background. It seems that a q — 8 
GeV 2 far 1 is enough to induce a dijet asymmetry comparable to what is observed in data and the differences between 
the curve with q — GeV 2 far 1 and q — 8 GeV 2 far 1 for the dijet azimuthal correlation are in agreement with 
the differences between the CMS Monte Carlo and CMS data. Moreover, Q-PYTHIA results show the same trend 




(a) Dijet Asymmetry. (b) Dijet Azimuthal Correlation. 

Figure 2: Jet Observables for a simulation using Q-PYTHIA with different q parameters, embedded in a background with T = 0.9 (cr ~ 1 1 GeV). 
The red/dotted curve corresponds to a q = GeV 2 fin , the black/dashed to a q = 4 GeV 2 fm - ' and the green/solid to a q = 8 GeV 2 fm -1 . The 
blue points are the CMS data while the purple triangles the CMS Monte Carlo simulation. The jet area method is used to subtract the background. 



as the average missing transverse momentum data measured by CMS[2 |, as one can see in figure [3] Our toy model 
simulates a system with global variables similar to those measured in the experiments: multiplicity, background level 
and background fluctuations. However, the distribution of particles in space and the spatial range of their correlations 
remains unconstrained. We believe that the missing p ± observable is extremely sensitive to this last feature, so in the 
absence of a truly realistic background model, we choose to present only the results for the Q-PYTHIA simulation 
without background while postponing more detailed studies which require truly realistic backgrounds. As reference 
we compare the simulations without quenching, like standard PYTHIA (q = GeV 2 far 1 ) and the CMS simulation 
for this observable. Our results for the missing transverse momentum outside a cone of AR = 0.8 around the leading 
or subleading jet axis, displayed i n figu re |3(a)| are in well agreement with the CMS simulation. By setting q — 8 
GeV 2 far 1 in Q-PYTHIA (figure 3(b) i we can compare our simulation with the CMS data. Again, we see that, at 
least qualitatively, it is able to describe the main features of the data, i.e., in general a softest composition of the event, 
and even outside the cone, the composition is still soft. It seems that a standard jet quenching model, where the 
energy loss is within a close range of the parent parton, is able to explain the soft composition outside the cone, as 
well as the other two observables, the dijet asymmetry and azimuthal correlation, without a compelling need of a large 
angle emission mechanism. This can be interpreted as a consequence of event topology since that already in standard 
PYTHIA, there are events with a large dijet asymmetry (Aj > 0.3), which are characterize by a large hard composition 
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(a) < p\ > with q = GeV 2 fin" 1 . 



(b) < p" > with q = 8 GeV fin- 



Figure 3: Average missing transverse momentum < p T > outside of a cone (AR > 0.8) around the leading or subleading jet axis as a function of 
the dijet asymmetry, A j for Q-PYTHIA simulation. Colored bands show the contribution to the < pjj, > for five ranges of particles pj. 



outside of the cone (see figure 3(a) I. This contribution can only come from an event in which the subleading jet emits 
a hard parton at large angle, originating itself a third semi-hard jet. When propagating through a medium, all jets, 
including this third jet, will be quenched, and consequently, their initial components are transformed into more soft 
momentum particles. 



5. Conclusions 

By comparing the two background subtraction techniques we saw that jet energy reconstruction shows sensitivity 
to fluctuations of the background when using FastJet. In turn, the CMS-like method has a dependency on the pa- 
rameter E T j ets and presents some deviations in the angular reconstruction of the dijet pair. This can be related to the 
intrinsic structure of the background, indicating that may be needed more than an effective p and <x to characterize 
a background. As for the description of the observables with Q-PYTHIA, it seems that this model goes in the same 
direction than data for the three observables addressed in this study: large jet momentum imbalance and an azimuthal 
correlation which does not deviate much from the proton-proton case for a reasonable value of q; and the presence of 
the higher amount of soft particles at large angle without the need of additional mechanisms of large angle emission. 
Acknowledgments: This work was supported by the European Research Council grant HotLHC ERC-2001-StG-279579; by 
Ministerio de Ciencia e Innovation of Spain grants FPA2008-01 177, FPA2009-06867-E and Consolider-Ingenio 2010 CPAN 
CSD2007-00042; by Xunta de Galicia grant PGIDIT10PXIB 206017PR ; by Fundacao a Ciencia e para a Tecnologia of Portugal 
under projects SFRH/BD/64543/2009 and CERN/FP/1 16379/2010 (LA); and by FEDER. 



References 

[1] G. Aad, et al., Observation of a Centrality-Dependent Dijet Asymmetry in Lead-Lead Collisions at sqrt(S(NN))= 2.76 TeV with the ATLAS 

Detector at the LHC, Phys.Rev.Lett. 105 (2010) 252303. |arXiv: lPlX . 6182 
[2] S. Chatrchyan, et al., Observation and studies of jet quenching in PbPb collisions at nucleon-nucleon center-of-mass energy = 2.76 TeV, 

Phys.Rev. C84 (2011) 024906. [arXiv : 1 102 ■ 1957] 
[3] S. Chatrchyan, et al., Jet momentum dependence of jet quenching in PbPb collisions at sqrt(sNN)=2.76 TeV, Phys.Lett. B712 (2012) 176-197. 

larXiv:T2 02.5022 

[4] N. Armesto, L. Cunqueiro, C. A. Salgado, Q-PYTHIA: A Medium-modified implementation of final state radiation, Eur.Phys.J. C63 (2009) 
679-690. larXlv: 0907. 10141 

[5] B. Abelev, et al., Measurement of Event Background Fluctuations for Charged Particle Jet Reconstruction in Pb-Pb collisions at y.Sjv/v = 2.76 

TeV, JHEP 1203 (2012) 053. |arXiv: 1201 .2423| 

[6] M. Cacciari, G. P. Salam, G. Soyez, FastJet user manual, Eur.Phys.J. C72 (2012) 1896. |arXiv: 1111 .6097| 

[7] O. Kodolova, I. Vardanian, A. Nikitenko, A. Oulianov, The performance of the jet identification and reconstruction in heavy ions collisions 
with CMS detector, Eur.Phys.J. C50 (2007) 1 17-123. 



